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We reveal by first-principles calculations that in iron pnictides Sr4M2 0eCrFeAs2 (M=Sc, Cr) the 
two-dimensional CrFeAs2 layers exhibit a robust band structure of half-metallic antiferromagnet 
(HMAFM). Due to the thick blocking layer, the interlayer coupling is vanishingly small and thus 
the conductive channels in individual layers may take the alternative spin direction randomly. We 
show that, since the spin magnetizations of Fe and Cr are different in a ferromagnetic state, applying 
a strong magnetic field and ramping it down gradually can align the spin direction of conductive 
channel of all HMAFM layers, which restores the bulk HMAFM. 

PACS numbers: 74.70.Xa, 85.75.-d, 31.15.A- 



Spin is one of the fundamental properties of elec- 
tron. Exploitation of the spin-dependent transportation, 
dubbed as spintronics, which emerged from the discov- 
ery of giant magnetoresistance effect [TJ [2] , has been a 
topic of particular interest in the last decades. Genera- 
tion of the spin-polarized current is the primary require- 
ment of spintronics devices. Half metals [3 j, a class of 
materials possess states in only one spin channel at the 
Fermi level (E-p), have received considerable attention 
since they yield fully spin-polarized currents. 

Owing to the asymmetric electron populations in the 
two spin channels, half metals are often ferromagnetic 
(FM) or ferrimagnetic. Half-metallic antiferromagnet 
(HMAFM) materials proposed first by van Leuken and de 
Groot[4 should be superior in many applications, where 
stray fields are harmful, since they achieve fully spin- 
polarized currents without showing net magnetization. 
To date, many HMAFM candidates have been engineered 
by the first-principles calculations based on the double 
perovskite compounds [5HZ] and Heusler alloys [8 . Be- 
sides, HMAFM have been predicted by doping in oxide 
semiconductor [9 , transition- metal (TM) oxides [TO] and 
cupratepj]. 

In a previous study [12 , the present authors proposed 
to achieve HMAFM based on the recently renaissant iron 
pnictides [TO] . We substituted half of the Fe atoms by Cr 
in BaFe2As2- The ground state was characterized by the 
following properties: the Fe and Cr atoms form inter- 
vening three-dimensional (3D) lattice, the spin magne- 
tizations of Fe, Cr and As compensate completely, and 
close to E F the states are contributed exclusively from 
the spin channel of Cr:3d majority electrons. The mate- 
rial BaFe2As2 is therefore predicted as a HMAFM. 

In the present work, we focus on the recently discov- 
ered iron pnicitides Sr 4 M 2 06Fe2As2 (M=Sc, Cr) [14]. 
These new iron pnictides (see also [TOJ [TO]) are unique 
in the structure where the FeAs layers are separated 
by layers of perovskite structure as thick as 9.8 A , 
and thus is much more two-dimensional (2D) than other 
iron pnictides. We substitute half of the Fe atoms by 
the conjugate element Cr as revealed in Ref.[12]. From 
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FIG. 1. (Color online) Lattice and spin orders of four de- 
generated configurations of TM layers in the ground state of 
Sr4Sc2 0eCrFeAs2. The As ions and blocking layers in the 
unit cell are omitted for a clear view. Configuration I and II, 
in which all Cr (and Fe) ions take the same spin direction, are 
predicted to be HMAFM, while III and IV are not. 



first-principles calculations it becomes clear that in the 
ground state the individual TM pnictide (CrFeAs2) lay- 
ers take the checkerboard order of Cr and Fe ions as- 
sociated with the antiferromagnetic (AFM) spin order, 
and exhibit HMAFM electronic structure. Due to the 
vanishing interlayer coupling, both lattice and spin dis- 
orders easily appear in the c direction, which suppresses 
the bulk HMAFM property. Interestingly, it is found 
that the spin magnetizations of Cr ions are much larger 
than Fe ions in the FM state. Therefore, a post-process 
of applying a strong magnetic field and ramping down 
gradually can restore the spin order along the c direction 
and thus achieve the 3D HMAFM. The present approach 
to achieve the bulk property based on the robust 2D one 
prepared in advance may shed light for future exploration 
of HMAFM. 

The calculations are performed using the projector 
augmented-wave (PAW) method [17] as implemented in 
the VASP simulation package p~8j [TO], with the gener- 
alized gradient approximation in PBE type [20] for the 
exchange-correlation functionals. The crystal structures 
of the parent materials revealed experimentally [14] are 
implemented in the present calculations for the doped 
system. Structure optimization is then performed us- 
ing 8x8x2 Monckhorst-Pack [21] k-point mesh until the 
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FIG. 2. (Color online) Total density of states (DOS) and 
those projected on the TM pnictide layer and blocking layer 
of Sr4Sc20eCrFeAs2, respectively, for the atomic and spin 
orders I and II in Fig. [I] 

TABLE I. Spin magnetizations of Cr, Fe and As ions and 
total spin magnetization (a) in states I and II shown in Fig.JT] 
and (b) in the FM state (units: /xb). The total spin magneti- 
zation is obtained by integrating charge densities of both spin 
channels. 



Conf. 


Cr 


Fe 


As 


M tot 


(a) 


2.83 


-2.68 


-0.08 


0.00 


(b) 


2.45 


0.56 


-0.08 


11.33 



stress of the supercell and the force of the ions are less 
than 1 kbar and 10" 3 eV/A. 

We begin with the identification of ground state of 
Sr 4 Sc2 06CrFeAs2 by calculating the total energies of 
possible configurations within the y/2 x y/2 x 2 super- 
cell. The checkerboard configuration of Cr and Fe ions 
associated with AFM spin magnetizations is of energy 
lower than other 2D configurations at least by 0.23 eV. 
The robust in-plane checkerboard order is consistent with 
the case of BaCrFeAs2 [12] , but in sharp contrast with 
the study for the '1111' system [22]. 

On the other hand, due to the diminished coupling 
across the thick blocking layer, the four configurations 
shown in Fig. [I] with different orders along the c axis 
possess almost the same energy (AE < 3 meV). 

Figure [2] shows the total density of states (DOS) and 
those projected to the block layer and CrFeAs2 layer, 
for the configurations I and II in Fig. [I] where the Cr 
(Fe) 3d majority electrons are spin- up (spin-down). The 
band near Ep (from -2 to 1 eV) is dominated by states 
from CrFeAs2 layers, and all the occupied states from the 
blocking layer are deep in the valence band. This result 
reveals the role of blocking layer as the charge reservoir. 
At Ey the DOS is exclusively contributed by the spin- 
up states, which consist of the Cr:3d majority and Fe:3d 




FIG. 3. (Color online) From (a) the state III in Fig. [T] (2D 
HMAFMs) to (b) a bulk HMAFM. Left column: lattice and 
spin orders; center column: DOSs of individual TM layers; 
right column: total DOS. DOSs are shown in energy range 
from -1 to 1 eV relative to Ef. 



minority electrons, and a part of As:4p states. As illus- 
trated in Table [I] the total magnetization is zero due to 
the complete compensation of the spin magnetizations of 
Cr, Fe and As ions. Therefore, the material is a HMAFM 
in the states I and II shown in Fig. [I] 

In BaCrFeAs2 [12] . the strong interlayer coupling via 
the direct As bonding delocalizes the As:4p states, result- 
ing a narrow band gap of 0.3 eV. In Sr4Sc20eCrFeAs2, 
the coupling is diminished due to the thick blocking 
layer. Therefore, the As:4p states are well confined in 
the CrFeAs2 layer. The energy gap of Sr4Sc20eCrFeAs2 
is calculated to be 0.7 eV, which is much larger than that 
ofBaCrFeAs 2 . 

For the configuration III in Fig. [T] the total DOS is 
symmetric in the two spin channels, as shown in the right 
column of Fig. [3ja). However, projecting the total DOS 
on the blocking layers and the two individual TM pnic- 
tide layers, we find that, first, the blocking layers still do 
not contribute to states close to Ef, and secondly, the 
two TM pnictide layers are HMAFMs with the metallic 
channels in opposite spin directions, as shown in the cen- 
ter column of Fig.[3ja). The same results are obtained for 
configuration IV in Fig. [I] Therefore, Sr4Sc20eCrFeAs2 
shows a robust 2D HMAFM, while the bulk HMAFM 
may be suppressed due to the degenerated states of dif- 
ferent spin orders along the c axis. 

Since the spin order can be tuned by external magnetic 
field, we investigate the property of FM spin configura- 
tion, which is achieved when the external magnetic field 
is strong enough. It is revealed by first-principles calcu- 
lations that in the FM state the spin magnetization of Cr 
is much larger than that of Fe as displayed in Table [TJ in 
a sharp contrast to the AFM state where they are almost 
the same. The difference in the spin magentizations, and 
thus in the Zeeman energies, can be used to tune the 
spin orders. In a strong magnetic field overcoming the 



3 




FIG. 4. (Color online) Band structure of Sr4Sc20eCrFeAs2. 
(a) spin-up band; (b) spin-down band, (c) Fermi surfaces of 
the spin-up band. 



energy difference between the FM and AFM states, the 
spin magnetizations of both Cr and Fe are simply aligned 
along the external field. During the process of ramping 
the external field down to zero, the spin magnetization 
of Fe flips to the direction opposite to the external field, 
since the spin magnetization of Fe presumes smaller Zee- 
man energy and the AFM state is much more stable than 
the FM state. Thus, no matter what the initial c-axis 
spin order is, the AFM order between Cr and Fe can be 
restored over the bulk material as shown in Fig. (3^b). 
The 3D HMAFM property is eventually achieved. 

We then look at the band structure of 3D HMAFM 
Sr4Sc2 0eCrFeAs2. For the configuration I in Fig. [T] 
we use a y/2 x y/2 x 1 supercell to perform the calcu- 
lation. As displayed in Fig. Qa), the partially occupied 
bands in the spin-up channel include two hole bands con- 
tributed by the TM:t2# electrons and two almost degener- 
ated Fe:d 3 ^2_ r 2 bands corresponding to the peak in DOS 
near Ep in Fig. [2] The Fermi surfaces of the spin-up 
channel are shown in Fig. [4^c). The hole pockets near 
the T point are induced by the TM:t2# states and two 
Fermi surfaces near the X points of the Brillouin zone are 
contributed by the Fe:d 3 ^2_ r 2 states. All the Fermi sur- 
faces show a more significant 2D character as compared 
with the iron-pnictides with thinner blocking layers [23] . 
There is no Fermi surface in the spin down channel since 
a gap is opened. 

First-principles calculations reveal that the band struc- 
ture of configuration II is almost the same as that of 




FIG. 5. Band structure of Sr4Cr2 0eCrFeAs2. (a) spin-up 
band; (b) spin-down band. 



configuration I, both characterized by the same Cr (Fe) 
spin direction over the bulk. Therefore, we predict that, 
as an aggregate of configurations I and II, the bulk 
Sr4Sc20eCrFeAs2 subjected to the processing of an ex- 
ternal field is of the same band structure as shown in 
Fig. [4] 

Next let us address whether Cr atoms in 
Sr4Sc206CrFeAs2 would go to the blocking layers. 
In the blocking layer, only Sc atoms might be sub- 
stituted by Cr, which forms the Cr^c-Sco anti-site 
defects. Using a 2^/2 x 2^/2 x 1 supercell we calculate the 
formation energy of anti-site defect, which is estimated 
to be as large as 2.41 eV. The Crs c -Sccy anti-site defect 
is thus effectively prohibited even the synthesizing 
temperature reaches 1500 K [14]. Therefore, the doped 
Cr ions prefer exclusively the Fe sites. 

Finally we turn to the twin parent material 
Sr4Cr20eFe2As2, where the Cr atoms exclusively occupy 
the B-sites of perovskite structure in the blocking layer 
[T4] . In order to manipulate the band structures, we 
try the same approach, namely replacing half of the Fe 
atoms with Cr atoms in the FeAs layers. We find that in 
CrFeAs2 sheets the Fe and Cr atoms prefer the checker- 
board distribution associated with an AFM order among 
Fe and Cr ions same as the case of Sr4Cr20eCrFeAs2, and 
that Cr ions in the blocking layers form AFM order them- 
selves consistently with the experimental observation 
[24] . As shown in Fig. [5j the material Sr 4 Cr 2 06CrFeAs2 
has a gap in the spin-down channel and is metallic in 
the spin-up channel, and thus is predicted to be another 
HMAFM. The gap is 0.2 eV, smaller than the previous 
case since Cr atoms in the blocking layers yield additional 
bands around E^. In the majority band, the orbitals d yz , 
d xz and d x 2_ y 2 of Cr atoms are fully occupied, while the 
other two orbitals are well above Ey due to the Coulomb 
repulsion between them and p electrons of oxygen; all 3d 
orbitals in the minority band are empty. 

To summarize, we reveal by first-principles calcu- 
lations that in iron pnictides Sr4Sc2 06CrFeAs2 and 
Sr4Cr206CrFeAs2 the CrFeAs2 layers exhibit a robust, 
two-dimensional property of half-metallic antiferromag- 
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net. Although the weak interlayer coupling may allow 
spin disorders along the c axis, the antiferromagnetic spin 
order between Cr and Fe can be restored over the bulk 
material by applying a pulse of external magnetic field, 
which gives rise to three-dimensional half-metallic anti- 



ferromagnet. 
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